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Abstract. We compare and contrast three different methods for the con- 
struction of the differential relations satisfied by the fundamental Abelian 
functions associated with an algebraic curve. We realize these Abelian func- 
tions as logarithmic derivatives of the associated sigma function. In two of 
the methods, the use of the tau function, expressed in terms of the sigma 
function, is central to the construction of differential relations between the 
Abelian functions. 



1. Introduction 

In the mid-1970s, the results of Its, Matveev, Dubrovin and Novikov (see 
[9]) led to the discovery of a remarkable 6'-functional formula to solve the KdV 
equation ut = 6uUx — Uxxx- This solution was given as a second logarithmic 
derivative of a Riemann theta-function: 

(1.1) u(x,t) = -—\n9(Ux + Vt + W) + C 

with U , V, VV = const G and C G C. The theta-function in that formula 
was constructed from a hyperelliptic curve Xg of genus g, while the "winding 
vectors" U , V are periods of abelian differentials of the second kind on Xg. 
Further, this formula is in a sense universal; it was generalized by Krichever 
[20] to other integrable hierarchies - whose solutions were associated with other 
algebraic curves. In this paper we consider the converse problem: 

Given an algebraic curve Xg, of genus g, its Riemann period matrix r, and 
its Jacobi variety Jac(Xg) = C/(lc,©r), we may construct 9-functions 9{u; t), 
u G Jac(Xg); then the fundamental Abelian functions on Jac(Xg) may be 
realized as the second logarithmic derivatives of9{u,T), pij{u) = — ^ du^du'^^ ■ 
We wish to construct all differential relations between these Abelian functions 
on Jac(Xg). 

The simplest case, the Weierstrass cubic, y'^ = Ax^ — g2X — g^, is an algebraic 
curve of genus one. This is uniformized by the Weierstrass elliptic functions. 
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X = p{u),y = p'{u), and the differential relations are 

(1.2) p" = 6p2-|, p'2^4p3-^2p-^73. 

In the case of higher genera, g > 1, the derivation of analogous equations be- 
comes much more complicated. In particular, the fundamental Abelian func- 
tions are now partial derivatives of a function of g variables. The different 
approaches to this problem form the main content of the paper. We restrict 
our analysis to the case of (n, s)-curves introduced and investigated in this 
context by Buchstaber, Leykin and Enolski [3l H] 



1.3) X,: y^ = x'+ A,,xy. 



ni+sj<ns 

These represent a natural generalization of elliptic curves to higher genera, 
and include the general hyperelliptic curve {n = 2). 

To any such curve we may associate an object which is fundamental to all 
our treatments of this problem, the fundamental bi-differential. This is the 
unique symmetric meromorphic 2-form on Xg x Xg, whose only second order 
pole lies on the diagonal Q = S, and which satisfies 

oo{Q,s) - ^^^Q^^_^^^^^^^^ = ms),aQm{Qms), 

where $(^(S'), ^(Q)) is holomorphic, and ^{Q),^{S) are local coordinates in the 
vicinity of a base point P, ^(-P) = 0. Usually co{Q, S) is realized as the second 
logarithmic derivative of the prime- form or theta-function [H]. But in our 
development we use an alternative representation of u{Q,S) in the algebraic 
form that goes back to Weierstrass and Klein, and which was well documented 
by Baker [I] 

(1.4) u{Q, S) = ,^,'^^,^;f ^ r^dxdz + 2du{Qf^du{S), 

fy{Q)fw{S}{x - ZY 

where Q = {x,y), S = {z, w), and the function J^{Q, S) = J^{{x, y), {z, w)) is a 
polynomial of its arguments with coefficients depending on the parameters of 
the curve Xg du = (dui, . . . , dug)'^ is the vector of basic holomorphic differen- 
tials. The factor 2 in the final term is introduced so that this theory reduces 
to the Weierstrass theory in the elliptic case. Here, x is a symmetric matrix 
expressed in terms of the first and second period matrices, 2u, 2// respectively, 
as X = u'^rj. This provides the normalization of uj{Q, S). We will refer to the 
first term on the right of (11.41) . which involves the polynomial J^{Q., S), as the 
algebraic part, so that 

(1.5) u{Q, S) = w^'*5(Q, S) + 2du{Q)^xdu{S). 
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This representation was revisited and developed by Buchstaber, Leykin and 
Enolski [3] and more recently by Nakayashiki [23] . 

The algebraic representation of the fundamental differential, as described 
above, lies behind the definition of the multivariate sigma function in terms of 
the theta-function. This differs from Q by an exponential factor and a modular 
factor: 

(1.6) a{u) = C(r)exp {^u^uj-^riu} 9 {^co'^u; r) . 

Here the g x g matrices 2u, 21] are the first and second period matrices, and 
r = u~^u'. The modular constant C(r) is known explicitly for hyperelliptic 
curves and a number of other cases. However its explicit form is not necessary 
here, for the fundamental Abelian functions are independent of C(r). These 
modifications make cT(ii) invariant with respect to the action of the symplectic 
group, so that for any 7 G Sp{2g, Z), we have: 

(1.7) a{u;jT)=a{u;T). 

The multivariate sigma-function is the natural generalization of the Weier- 
strass sigma function to algebraic curves of higher genera, i.e. {n, s)-curves in 
this context. In his lectures [30], Weierstrass started by defining the sigma- 
function in terms of series with coefficients given recursively, which was the key 
point of the Weierstrass theory of elliptic functions. A generalization of this 
result to the genus two curve was started by Baker |2j and recently completed 
by Buchstaber and Leykin [6j, who obtained recurrence relations between coef- 
ficients of the sigma-series in closed form. In addition, Buchstaber and Leykin 
recently found an operator algebra that annihilates the sigma-function of a 
higher genera (n, s)-curve [7j. The recursive definition of the higher genera 
sigma-f unctions remains a challenging problem to solve, with [7J providing a 
definite step. We believe that the future theory of the sigma and corresponding 
Abelian functions can be formulated on the basis of sigma expansions that will 
complete the extension of the Weierstrass theory to curves of higher genera. 

In this paper we study the interrelation of the multivariate sigma and Sato 
tau functions. The r-function was introduced by Sato [281 129] in the much 
more general context of integrable hierarchies. But it seems there are few re- 
sults studying algebraic curves in Sato theory. However we should mention 
recent work by Konopelchenko and Ortenzi analyzing algebro-geometric struc- 
ture in Birkhoff strata of the Sato Grassmannian [19J. The recent papers of 
Matsutani and Previato [211 [22], studying Jacobi inversion on Jacobian strata 
of (r, s) curves, is also relevant to our work, relating stratification of the Sato 
Grassmannian to partitions. 

Here we deal with the 'algebro-geometric r-function' (AGT) associated with 
an algebraic curve. 
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The AGT of the genus g curve Xg is defined, following Fay, [I5l HE], as a 
function of the 'times' t = {ti, . . . ,tg,tg+i, . . .), a point u G Jac(Xg), as well 
point P G X^; it is given by the formula 

r(t; U,P) = Uk{P)tk + ^^'^'^ ) 1 ^ E ^mn{P)tmtn 

\k=l J I m,n>l 

Here the "winding vectors" Uk{P) appear in the expansion of the normalized 
holomorphic integral v, the quantities Umn{P) define the holomorphic part of 
the expansion of the fundamental differential of the second kind u{Q, S) near 
the point P. We then introduce the r-function by the formula: 

This representation of r in terms of a was used by Harnad and Enolski [13j 
to analyze the Schur function expansion of r for the case of algebraic curves. 
Recently Nakayashiki [25] has independently suggested a similar expression for 
the AGT in terms of multivariate a-functions. In this paper we concentrate on 
the application of this representation to the derivation of the differential rela- 
tions between Abelian functions of the (n, s)-curve, continuing and developing 
the work of |13j . 

Developing a further analogy with the Weierstrass theory of elliptic func- 
tions, we represent the Abelian functions, that is, 25f-periodic functions 

F{u + 2nuj + 2nV) = F{u), Vn, n' G N 

on the Jacobian 

Jac(Xj = C/2uj © 2uj' = ^"^Jac(Xg), 2u = A, 
as second and higher logarithmic derivatives: 

(1.9) au) = ^\na{u), 

Q2 Q3 

(1.10) pij{u) = --^—^\ncr{u), pijk{u) = .g^^ ^^ ^('^)' ^^^-^ 

where i,j,k,--- = 1 . . . ,g. We should remark that the Qi'^) S'-i's not Abelian 
functions. In this notation, the genus one Weierstrass equations (11. 2p become 

pun = Qpli - y, Piii^ = - 5-2^11 - 93- 

For general g, the pij, pijk, ■ ■ ■ are called Kleinian p-functions. They are con- 
venient coordinates to represent the dependent variables in the hierarchy of 
integrable systems. 
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In this paper we compare and contrast three approaches to obtain the partial 
differential relations for the Abelian functions associated with the (ra, s)-curve 
Xg, using Kleinian p- functions as coordinates. 

The first of these, and the best known, is the classical approach of comparing 
two different expansions of the fundamental bi-differential; this yields first the 
solution of the Jacobi inversion problem for the curve, and in higher orders, a 
sequence of differential relations involving the pij. 

The r-function approach to the derivation of completely integrable systems 
of KP type has led to two different ways [8J to obtain relations between Taylor 
coefficients of the r function expansion. The first of these specifically exploits 
the fact that these Taylor coefficients are determinants, i.e. Pliicker coordinates 
in the Grassmanian, and they hence satisfy the Pliicker relations [15]. The 
second is based on the Bilinear Identity, which leads to the Residue Formula 
[T6j , giving a family of differential polynomials in r which must vanish - these 
are the partial differential equations we need. 

We further consider and compare the two techniques based on the r-function 
method, which give a derivation of the required differential relations; spe- 
cializing to a particular algebraic curve, we consider its algebro-geometric r- 
function. The special feature of our development is that we define this r- 
function in terms of the multidimensional a-function of the curve, leading 
to coordinates that are explicitly written in terms of Kleinian p-functions. 
The differential relations we find between these functions can be understood 
as arising from special solutions of integrable hierarchies of KP type, associ- 
ated with the given curve (see for example [H [3l [23]). We will describe the 
correspondence between individual differential equations for p-functions with 
Young diagrams defining Pliicker relations. We illustrate these approaches by 
considering two particular examples: the genus 2 hyperelliptic curves and 
the genus 3 trigonal curve, which can be found in [21 El fT^). 

We remark that in similar context, Nimmo [26] in 1982 applied the methods 
of symmetric function theory to describe the KP hierarchy, and studied the 
action of the recursion operator. The general approach based on Pliicker co- 
ordinates for deriving KP-flows in terms of Kleinian a-functions was recently 
discussed in [13]. Here we develop these ideas and consider some non-trivial 
examples to clarify the interrelation of the r-functional formulation of inte- 
grable hierarchies and the cr-functional approach. We will also consider the 
relationship between of this derivation, based on the Pliicker relations, and 
that based on the Residue Formula. Both give a systematic way of generat- 
ing the required relations, but the differences between the two approaches are 
instructive. 
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2. Algebraic curves 
Let Xg be a genus g > 1 algebraic curve given by the polynomial equation 
(2.1) fix, y) = 0, fix, y) = y'' + y"-^ai(x) + . . . + ao(x). 



We shall consider in what follows two relatively simple curves of the class (12. ip . 
Example I: the hyperelhptic genus two curve 



(2.2) 
and 



y = 4x + a/^x + . . . + ao 



Example II: the cyclic trigonal genus three curve 



(2.3) 



fix, y)=y^ - (a;^ + /igx^ + /igX^ + /igX + /i 



12 



We equip Xg with a canonical basis of cycles (oi, . . . , a^; bi, . . . , b^) G HiiX, Z) 
We denote by du = (dwi, . . . , dug)^ the vector whose entries are independent 
holomorphic differentials of the curve Xg as well as their a and b-periods, 



(2.4) 



2uj 



dui 



2uj' 



dui 



The period matrix (20;, 2u') is the first period matrix, and the matrix r 
u~^u' belongs to the upper Siegel half-space, (3 : = r, Imr > 0. 



The 6'-function 6[a]iz;T) with characteristics [a] 



a 



, [2a] eZ3 xZ3 



of the algebraic curve Xg of genus g is defined through its Fourier series 
e fr (z-r) 



exp |m (m + a)'^ r (m + a) + 2m (m + a)'^ (z + 6) | 



and possesses the periodicity property: for arbitrary a,b E and arbitrary 
a', b' G the following formula is valid 
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z + ra' + b'; r) 

--e[ia + a' fib + b'f ] iz; r)exp |-ma'^ra' - 2ma'^z - 2m(b + b')^a'} 
We introduce the associated meromorphic differentials dr = (dri, . . . , dvgf 



and their periods 



2r] 



dvi 



2rl 



dri 
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which form the second period matrix (2?7, 1vl\ The period matrices satisfy the 
condition 



?7 T]' J \ — Ig J \ V V' J 2 \ — Ig 

Here we denote the half-periods of the holomorphic and meromorphic dif- 
ferentials by (a;, u') and [r], rj') in order to emphasize the analogy with the 
Weierstrass theory. We will also use the notation A = 2u and B = 2u' for the 
periods of holomorphic differentials. Further we denote 

(2.5) dv{Q) = (dviiQ), . . . , dvg{Q)f = A-^du{Q) 

as the vector of normalized holomorphic differentials. 

The explicit calculation of canonical holomorphic differentials and the mero- 
morphic differentials conjugate to them is well understood; in particular we 
have 

Example I 

xdx dx 
dui = , du2 = — , 

y y 

x"^ dx x(a3 + 2aiX + 12x^) dx 
dri = , dr2 = 



Ay 



Example II 



, dx - xdx , dx 
dui = — , duo = — TT, du-i = — -, 
Sy iy^ iy^ 

x^ dx 2xdx (5x^ + 3;U3X + yUg) 
dri = — dr2 = — , drg = dx. 

Remark 2.1. Note that our labelling of the differentials is the reverse of [T2] . 
with the interchange 1 ^ 2 in example I, and (1,2,3) ^ (3,2,1) in example II. 

We introduce the fundamental bi-differential uj{Q, S) on X x X which is 
uniquely defined by the conditions: 

(i) it is symmetric: 

(2.6) Lo{Q,S) = tuiS,Qy, 

(ii) it has its only pole along the diagonal Q = S, in which neighbourhood 
it is expanded in a power series according to 

(2-^) ^) - (Hrfil)^ ^ S ^™n(P)e(Q)"-^e(5)"-Me(g)de(S); 
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(iii) it is normalized such that: 
(2.8) {uj{Q,S) = 0, j = l,...,g. 



The well known realization of the differential uj{Q,S) involves the Schottky- 
Klein prime form E{Q,S), which is a (— 1/2, — l/2)-differential defined for 
arbitrary points Q, S & X 

where is a 6'-function with non-singular odd characteristics [a], dv is 

the vector of normalized holomorphic differentials, and 

a o 

k=l ^ 

The bi-differential ^(Q, S) is then given by [H]. 

(2.10) uj{Q,S) = dQds\YiE{Q,S). 

We emphasize that in this paper, we will instead rely on alternative "alge- 
braic" constructions of the differential u{Q,S). By following classical works 
such as [milH], together with results documented in [1] we express the differ- 
ential uj{Q, S) in the form 

(2.11) coiQ, S) = r^dxdz + 2 du{Q)^ xdu{S) , 

fy{Q)fw{S){x - zy 

where Q = {x,y), S = {z, w), and the function J^{Q, S) = J^{{x, y), {z, w)) is a 
polynomial of its arguments, with coefficients depending on the moduli of the 
curve Xg. Finally, x is a symmetric matrix x'^ = x that is chosen to provide 
a normalization of u{Q,S); it is expressible in terms of the first and second 
period matrices x = uo'^rj. We will refer to the term of uj{Q, S) including the 
polynomial J^{Q, S) as its algebraic part, 

(2.12) w(g, S) = w^^'^(g, S) + 2 du{Qfxdu{S) 

In the vicinity of a base point P, where points Q and 5* are represented by local 
coordinates ^(Q) and ^{S) respectively, the holomorphic part of ijJ^^^{Q,S) is 
expanded in the series 



J-'{Q, S)dzdx 



daQ)das) 



fy{Q)fw{S){x - Z 

(2.13) 

= E<f(^)^(^)'^(^)'de(g)de(^). 

k,l=0 
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An algorithm to construct the polynomial J^{Q, S) is known, see e.g. [1] 
and therefore functions such as ujff{P) are needed for our construction are 
considered as known. In all that follows we will take the fixed base point P to 
be (oo, oo). We shall present below some explicit expressions for J-' as well as 
the first few terms of the expansions uff{P) in the simplest cases: 

Example I: 

(2.14) uj^{Q,S) = — , Q = {x,y),S = {z,w), 

A[x — z)'^ y w 

where 

(2.15) F{x,z) = 4:x'^z'^{x+z)+2a4x'^z'^+a3Xz{x+z)+2a2XZ+ai{x+z)+2ao- 

The polynomial F{x,z), with the properties F{x,z) = F{z,x) and F{x,x) = 
2y'^, is sometimes called the "Kleinian polar". 
Expanding this about the base point gives: 

ale "4 



0,0 - g , 
, alg _ , alg _ r, 



alg _ , alg _ 16a3 — 

128 

. . . etc. 



^0,2 — ^2,0 



I ,alg _ n 

^1,1 — u, 



Remark 2.2. For hyperelliptic curves, the coefficients uf^ vanish if either of 
i or j is odd. 

Example II: 



(2.16) 



u^Hi^,y),{z,w)) 



J-'{{x, y), {z, w)) dxdz 



{x - zyfy{x,y)f^{z,w) 
with the polynomial J-'((x,y); {z,w)) given by the formula 
J^{{x,y), {z,w)) = w^y^ 



(2.17) 



+ w \ w 



f{x,y) 

y 



+ T{x,z) ] + y iy 



w 



and 

(2.18) T(x, z) = 3/ii2 + {z + 2x)fig + x{x + 2 z)hq 

(2.19) + S/Xgx^^ + x^z^ + 2 x^z. 
Expanding this about the base point gives: 



,alg 
0,0 



0, 



T(z,x] 
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alg _ , alg 



2 



alg alg 2 I 5 2 

alg _ , alg _ 2 ,, I ^ ,,2 



. . . etc. 

Remark 2.3. uf^ = unless i + j = ImodS. This is a consequence of the 
cyclic symmetry of the curve. 



3. Algebro-geometric 6, a and t-functions 

Let V G Jac{Xg) and 6{v) be a canonical 6'-function, that is, a 6'-function 
with zero characteristics: 

^^(i;) = exp [inm^Tm + 2z7rt;^m} . 

The starting point of this paper is the transition from 6 to a-functions. For 
any point u G Jac(X) we define: 

(3.1) a{u) = C{T)9{A'^u)exp j^w^x-uj , 

where 6{v) is the canonical 6'-function, and C(r) is a certain modular constant 
that we do not need for the results that follow. We note that this a-function 
differs from the "fundamental a function" of the publications mentioned in 
the introduction by the absence of a shift of the 6'-argument by the vector 
of Riemann constants. Thus cr(0) 7^ 0. The sigma-function inherits a quasi- 
periodicity property from the corresponding 6'-function. 

We introduce the Kleinian multi-variable ( and p-functions as above (11.91 
ll.lUp . These functions are suitable coordinates to describe Abelian functions 
and the KP-type hierarchies of differential relations between them. 

In higher genera, the relations between Abelian functions become somewhat 
lengthy. These relations can often be summarized concisely by developing a 
matrix formulation of the theory, as was done |3] in the hyperelliptic case. 

The Sato-Fay algebro-geometric r-function of the genus g curve Xg of argu- 
ments t = {ti, ■ ■ ■ ,tg, tg+i, . . .)^, u G Jac(X), P G is defined as 

(3.2) r{t- u,P) = 9iY. Uk{P)tk + A-'u exp <^ - ^ 

^mn (-P) tmtn 

\fc=l / ^, m,n>l , 



a- AND r-FUNCTIONS 



11 



Here A is the matrix of periods of canonical holomorphic differentials and 
the winding vectors defined in (13.31) . Uk{P), appear in the expansion of the 
normalized holomorphic integral v in the vicinity of the given point P G X, 



dviQ')= dviQ') + J2Uk{P)aQ)' 

with ^ (Q) being the local coordinate of the point Q in the vicinity of the given 
point P, so that ^(-P) = 0. The quantities Umn{P) define the holomorphic part 
of the expansion of the fundamental second kind differential u}{Q, S) near the 
point P according to (12. 7p . 

We restrict ourselves to the case of algebraic curves with a branch point at 
infinity, and take P = (oo, oo) to be the base point where we expand all our 
functions. The winding vectors are in this case 

(3.3) Uk{oo) = Uk = A-'Rk, k = l,...,g, 

where -Ri, R2, . . . are residues of canonical holomorphic integrals multiplied by 
differentials of the second kind with poles of order k at infinity, giving: 



1 d'^-i 1 



Q= 



k d^Q)'-' '^'daQ) 

Using the above definitions, we can see that the algebro-geometric r-function 
is given by the formula, equivalent to ( 11. 8^ : 

(3.4) ^ = - — i;::- ' - exp ^ ^ > : 

In ( 13. 4p . uff is the algebraic part of the holomorphic part cu^ / of the expansion 
of the bi-differential as defined in ( I2.13p . One can see that the non-algebraic 
part, i.e. the normalizing bi-linear form, is absorbed into the a-function. 

Once we have an expression involving derivatives of the sigma function, we 
need to convert this to an expression involving derivatives of the p-function. 
To do this we start with the definition of the ( function ( II. 9p . which we write 
in the form 

ai{u) = Ciiu)a{u), i = l,...,g 

then repeated differentiation gives us a ladder of relations which enable us to 
recursively express any derivative of a in terms of pij,,,k, Q all evaluated 

at u. 

(3.5) Oijiu) = aj{u)Ci(u) - a{u)pij{u) 

cTijkiu) = ajk{u)Ci{u) - (y{u)pijk{u) - (Tk{u)pij{u) - (jj{u)pik{u), 
etc . . . 
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In the following sections, we will look at different methods for constructing 
such relations between the derivatives of a, and hence between the Abelian 
functions associated with the curve. 

4. The "classical" method 

The starting point for this approach is the Klein formula, which compares 
two different expressions for the fundamental bi-differential: 

Theorem 4.1. Let the canonical holomorphic differentials of the curve f{x, y) = 
be represented in the form 

Ui{x,y) , 
dui{x,y) = — -dx, i = l,...,g, 

fy{x,y) 

where Uk{x, y) are monomials of their variables. Then the following g relations 
hold 

k = l,...,g. 



(4.1) *J=i " " 



k=l 

J^{x,y; Xk,yk) 



{X - 

with polynomial J-'{P, Q) = J-'{x, y; z, w) defined in Ii2.12\) . an arbitrary base 
point of the Abel map Pq, and the corresponding vector of Riemann constants, 

As an application of this theorem consider the case of a curve with branch 
point at Po = (oo, oo). Let us tend P = (x, y) — )■ (oo, oo). Both, right and left 
had sides of (14. ip have poles at infinity. Equating principal part of poles, we 
get a set of relations between P^ = {xk,yk) and the multi-indexed p-symbols 
which can be interpreted as differential equations for the pjt,rfunctions. In this 
way, the solution of the Jacobi inversion problem can also be derived in terms 
of p-f unctions. 

4.1. Example: Hyperelliptic curve of genus two. The cr-functional real- 
ization of hyperelliptic functions of a genus two curve has already been dis- 
cussed in many places, see e.g. [H |3]. But we shall briefly describe here the 
principal points of the construction to convey the structure of the theory that 
we wish to develop for higher genera curves. 

We consider the genus two hyperelliptic curve of example I. The algebraic 
part of the fundamental bi-differential is given as above in eqn. (I2.14p [1]. 

Following the procedure in those papers, we equate principal parts of the 
highest (second order) pole to obtain 

(4.2) pi2 + xpii-x^ = 0, 
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which is the x-part of the Jacobi inversion problem for this curve. Then equat- 
ing principal parts of the next, lower pole and using (14. 2p leads to the relation 

(4.3) = -pii2 - 2;fcpiii, 

which completes the solution of the Jacobi inversion problem. 

We can use these relations to eliminate y and quadratic and higher terms in 
x, and we get from the next term 

This equation holds for both x^, so the coefficients of different powers of 
must be identically zero, and solving for pmi and piii2 we find 

(4.4) pun = 6pii + a4pii + 4pi2 + |a3, 

(4.5) piii2 = 6piipi2 + Q;4pi2 - 2p22, 

the first two 4-index relations in the genus 2 case. A relation for 3-index 
functions pjjfc can be obtained by cross-derivation of (14. 4p and (14.51) : 

Pl22 + P11P112 - P12P1II = 0. 

Other equations are found from higher order terms. At every stage we need 
to substitute for higher derivative terms (such as pnm, for example) by using 
derivatives of previously derived relations. In addition, multiplication by a 
3-index pij^ is sometimes useful, followed by substitution of known relations 
which are quadratic in the pijk- The first such relation is found to be 

Jace : pin = ^pX^ + "spii + a4pn + 4pi2pii + ^2 + 4p22- 

This is one of the relations that describe the Jacobi variety of the curve as 
algebraic variety. These relations are well documented in [2], see also [3]. 
Below we show an alternative way to derive these relations. 

We remark briefiy that one variation on the classical method is to invoke 
the use of the sigma function expansion in the Ui variables. The first few terms 
in the sigma expansion are derived from the lowest order expansion terms of 
(14. ip . These can then be used in a bootstrap fashion to derive higher order 
PDEs for the p functions. This approach was essential in the genus six (4,5) 
case considered in [lOj, which is discussed in more detail in that paper. 

5. Derivation of integrable hierarchies via Plucker relations 

The key to this approach is Sato's formula. Theorem (5.1) below. This gives 
an expansion of a ratio of r-functions, r(t; M)/r(0; u), as a series of rational 
expressions in the r-function and its derivatives att = 0. The i-dependence is 
given in terms of Schur polynomials s\{i) in the times ti. The coefficients of 
these polynomials are determinants of differential expressions of r, which are 
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Pliicker coordinates on a Grassmannian. Such coordinates satisfy the Pliicker 
relations - each partition A can be expanded in hooks, and the corresponding 
Pliicker coordinates are expressible, analogously to Giambelli's formula, as 
determinants of single hook partitions. These relations give the differential 
relations for the Abelian functions which we seek. 

For any partition A : ai > 0:2 of |A| = ^"=1 a^, the Schur 

polynomial of n variables xi, . . . , is defined by 

s^{x) = det (Pa,-i+i(a^))ij=i,„„„ , 

where the elementary Schur functions Pm{x) are generated by the series, 

00 f 00 

Pm{x)t"' = exp <^ ^ Xnf 
m=0 y n=l 

The first few Schur polynomials are 

Si{x) = Xi, 

S2{x) = X2 + \x\, = -X2 + \xl, 

S^{x) = X'i + X1X2 + \xl, S2^l{x) = -X'i + \xl, Si^i^i{x) = X2, - X1X2 + \xl 

Si{x) = 0:4 + X1X3 + \x\ + \x\x2 + ^x\, etc. 
The Cauchy-Littlewood formula 

exp <^ YnXnVn \ = Ysx{x)sxiy), 

ln=l j A 

where Sx{x) is the Schur function of the partition A : ai > 02 > . . . > 
leads to the Taylor expansion in the form 



f{x) = exp 




y=o 



y=0 A 

We now introduce the Frobenius notation for partitions. Single hook par- 
titions are denoted (a + 1, 1^). All partitions can be decomposed into 
finitely many hooks: 

(A) = (ai, . . . . . . ,/3r)- 

The total number of hooks r is called the rank of the partition A, [27J. In 
particular, a partition which decomposes into two hooks (rank 2) is written as: 

(r^,m,2^l') = {n - l,m - 2\k + I + l,k), 

where n>'m>l, k>0 and / > 0. Most of the formulae below are derived 
using rank 2 partitions. 
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Giambelli's formula now shows how to expand a Schur function s\ in hooks, 
sx{x) = det 

Theorem 5.1 (Sato formula ). Let T{t;u) be any function of vector argu- 
ments 

where Ui, U2 ■ ■ ■ is an infinite set of constant complex vectors from and 
u E is a parameter. Suppose that r(0; n) 7^ 0. Then for any partition 
(A) = («!, . . . , ar\/3i, . . . , f3r) and any u E 

(5-1) = E^AWdet , 

where the A(^rn\n){u) with m > 0, n > form a linear basis of the Grassman- 
nian: 

^(■m\n){u) = — yl(„|m)(— n) 

= (-ir+is„+i,in(5t)r(t;n)|^^Qr(0;n)-i 



(5.2) 



and 



^Pn+a+i{-dt)p {dt)T{t;u) 



a=0 



r(0;n) 



t=o 



• ,«r|/3l, • • • 



r_d_ l_d_ l_d_ 
* ~ [dfi' 2dr2' Sdfs' " ' ' 

Theorem 5.2 (Pliicker relations). For any partition (A) = (ai 
and any u G C^, the r-function satisfies: 

(5.3) TiO;uY-\s,idt)T{t;u)l^o = ^)L=o) • 

In particular, for any curve X, its Sato algebro-geometric r-function has the 
expansion 

rit: u) 



(5.4) 



r(0;n) 



= 1 + A{o\o){u)si{t) + A(i|o)(n)s2(t) 
+ A(o|i)(n)si,i(t) + . . . 

= 1 + A(o|o)ti + A(i|o)(i2 + |t?) + A(o|i)(n)(-t2 + |ti) + . . . 



where A(m|n) (1*) as defined in (15. 2p are the elements of the semi-infinite matrix 
A, 



A 



( ^(0|0)('") ^(0|l)(^i) ^(0|2)(li) 

A(i|o)(n) A(i|i)(w) A(i|2)(n) 

v4(2|0)(n) A(2|l)(w) ^(2|2)(W) 



(^0,^1, 
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with infinite vectors Aj, i = 1,2,.... To more complicated partitions such 
as (mi, . . . , mfe|ni, . . . , n^) we associate according to the Giambelli formula 
certain minors of A, 

^(mi,...,mfc|ni,...,nfc)(w) = det {^{mi\nj){u)) ^j^^ ,^ ■ 

Note that these A's are not independent - they satisfy the Pliicker relations. 
These are a family of differential equations satisfied by r, that represent a 
completely integrable hierarchy of KP-type. 

The definitions and relations given above are valid for any multivariate func- 
tion T{t; u). Below we consider functions r-functions constructed on the Jacobi 
varieties of algebraic curves. 

To each symbol we shall put in correspondence its weight 

9 

P^,2_2,...,g_3 ^ ^kjWj, 

fc-]^ ^2 kg 3 

where Wi is the order of vanishing of the holomorphic integral J dui at infinity. 
In other words, if the curve has a Wcicrstrass point at infinity, then Wi = 1 < 
W2 <■■■< Wg is the Weierstrass gap sequence at infinity. 

For a given weight W consider all Young diagrams which decompose into 
only two hooks, and write corresponding relations between multi-index symbols 
pii,....in- The first non-trivial Young diagram corresponds to the partition A = 
(2,2). There is only one multi-index symbol of weight 4, that is piiii{u) and 
its corresponding Pliicker relation for any curve is of the form 

(5.5) piiii{u) = polynomial of even symbols pij{u) 

For higher weights a larger number of multi-index functions can be constructed, 
and this number grows rapidly with increasing weight. To find them in terms in 
the form of polynomials of two-index functions, we shall write Pliicker relations 
corresponding to independent diagrams of the same weight, and solve the 
corresponding linear systems. The technique is best illustrated by examples. 

We consider the genus two hyperelliptic curve, with the differentials and the 
bi-differential chosen as in example I above. The simplest class of non-trivial 
Pliicker relations is found from the class of Young diagrams which may be 
decomposed into two hooks, that is, those of the form (2 + m, 2 + n, 2^, l'), 
with m > n > 0, A; > 0, / > 0. In Frobenius 'hook' notation these are (m -|- 
l,n\k + I + 1,1). The corresponding Pliicker relations read 

(5.6) 

These equations are all bilinear partial differential equations in the r-function. 
They may be expanded in terms of the Kleinian p^j and Q functions. 



A(n,k+l+l) ^{n,l) 
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The first Young diagram leading to a non-trivial Pliicker relation corresponds 
to the partition A = (2, 2), with Young diagram 



Writing (I5.3p in this case we obtain after simplification 

(5.7) KdV4 : piiii(u) = 6pii(u) + 4pi2(w) + a4pii(u) + ^03. 

This is of course the same equation as fl4.4p . The weight of the diagram is 
4 in this case, which is the same as the weight of the equation, defined as 
the weighted sum of indices in which the index "1" has weight 1 and index 
"2" has weight 3. We will use this weight correspondence in other cases too, 
and will denote the weight of the object by subscript i + 3j where i and j 
are respectively the numbers of I's and 2's in the multi- index relation. The 
relation fl5.7p is the analogue of Weierstrass' equation for in the genus 1 
case. 

The next group of diagrams are of weight 5, and correspond to the parti- 
tions A = (3, 2) and A = (2, 2, 1), with their transposes, which give the same 
equations, in the hyperelliptic case. Both these Pliicker relations lead to the 
equation 



+ 



KdV. = 



i.e. the Young diagram of weight 5 gives no new equations and we conclude 
that the following correspondence is valid 

A = (2, 2) 
A = (3,3), A = (2, 2,1) 

At weight 6 we have three independent Young diagrams with (2, 2) centres 



KdV. 



and 



The other two diagrams of weight 6 



and 



again give the same results as their transposes. 

Remark 5.3. The Pliicker relations associated with any Young diagram and 
its transpose are always the same for a hyperelliptic curve; however for general 
curves, this is no longer true. 
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These three independent diagrams give an overdetermined system of three 
equations. After substituting for higher derivatives of the known relation (15.71) . 
we can solve for the two unknowns puu and pf^ to get 

(5.8) KdVe : pui2 = Qpupu - 2p22 + a4pi2 

(5.9) Jace : plu = "^Pii + (^spu + o^aPu + 4:pi2pii + a2 + 4:p22 

The relation fl5.8p is another generalization of the equation for p" in the genus 
one case, and is identical to (14. 5 p found using the classical method. The 
relation (15. 9p corresponds to Weierstrass' equation for {p'Y in the genus one 
case. In what follows we will always assume that higher derivatives of the 
known relations at a lower weight have been eliminated. 

At weight 7 we have 5 independent Young diagrams with (2,2) centres. The 
resulting overdetermined system of equations, of weight 7, contain d multiplied 
by linear combinations of the two weight 6 relations (I5.9|5.8p . In addition we 
get a relation at weight 7, namely 

Puplll - PV12 - P\V2p\\ = 0, 

as noted above. We refer to relations of this type, linear in the three-index 
pijk, as quasilinear. These have no counterpart in the genus 1 theory. They 
can also be derived by cross-differentiation between the two relations (I5.7|5.8p . 
since 

d _ d 

— pun — ^ — piii2- 

OU2 OUi 

Equations KdVi and KdV^ describe the genus 2 solutions of the KdV hierar- 
chy associated with the given curve. To describe the Jacobi variety and the 
Kummer variety we must consider diagrams of higher weights. 

For the diagram of weight 8 we find a set of eight overdetermined equations 
with solution given by 

Jacg : p\^^ = oiQ- Ap22pi2 + a4p?2 + 4piiPi2 
KdVg : P1122 = 2piip22 + 4pi2 + |«3pi2, 
And at weight 9 the only new relation is the quasilinear relation 

8pl22pll - 4pi2pll2 - 4p222 + 2^4^122 - 03^112 " 4piiip22 = 

At weight 10 we have 18 overdetermined system of equations, which can be 
solved for the 3 functions of weight 10, P1222, P1121 P111P122, giving 

KdVio : P1222 = 6pi2p22 + "2^12 - - "o 

Jacio : P111P122 = + 2p22p\i + ^pupl^ + 0:2^12 + 4p22pi2 + l^spupii 

(2) 2 2 2 

Jacio : P112 = ao - 4p22pi2 + a4pi2 + 4piipi2 
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The equations Jacg, Jac^o , Jac^g^ represent an embedding of the Jacobi variety 
as a 3- dimensional algebraic variety into the complex space whose coordi- 
nates are pn, pi2, P22, piii, P112 

At weight 12 we get the final 4-index relation for P2222 and two quadratic 
3- index relations for P112P122 and ^111^222- We can continue in this manner 
at weight 14 to get more quadratic 3-index relations [3J. At the odd weights 
11, 13, 15, we get quasilinear relations which can also be found by cross- 
differentiation. As a practical point we note that the equations we derive can 
often contain ideals generated by the lower weight relations, and some work is 
required to identify genuinely new relations. 

At weight 16 we have 117 independent relations giving an overdetermined 
system of equations (we have checked only a selection of these). At this weight 
a new feature occurs. As well as the equations expressing the quadratic 3- 
index term P122P222 in terms of cubics in the pjj, we have terms which are 
quartic in the pi-,-. We can pick one of these quartic terms, say ^^2; ^^cl solve 
for this and for P122P222 to give us two relations. The relation involving a 
quartic in the is just the Kummer variety of the curve. This is the quotient 
of the Jacobi variety, Kum(X) = Jac(X)/(tt — )■ —u). In the case g = 2, the 
Kummer variety is a surface in C'^ which is given analytically by a quartic 
equation. This relation can can also be found from the identity 

(5.10) (Pm)(Pn2)-(PiiiPii2)' = 0. 

The same quartic also appears, multiplied by various factors, at higher weights. 

5.1. Example: Trigonal curve of genus three. As before we consider 
Example II, (2.3) whose holomorphic differentials are given in (2.x) [TT| I5| 
Here we generally follow the notation of [12j). The fundamental second kind 
differential is f ETTj) . 

The first Young diagram leading to a non-trivial Pliicker relation, as in the 
genus 2 case, corresponds to the partition A = (2, 2). Writing (15.31) in this case 
we obtain, after simplification, 

Pun = 6pii - 3p22- 

The weight of the diagram is again 4 in this case, which is the same as the 
weight of the equation, defined as the weighted sum of indices in which the 
index "1" has weight 1, index "2" has weight 2, and index "3" has weight 5. 
We will use this weight correspondence in other cases too, and will denote the 
weight of the object by subscript i + 2] + 5k where i,j and k are respectively 
the numbers of I's, 2's and 3's in the multi-index relation. 

In the trigonal case we no longer longer have the symmetry about the di- 
agonal of the diagram that we have in the genus 2 case, but we can restrict 
ourselves by taking the symmetric or antisymmetric combination of the two 
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diagrams related by transposition. In the weight 5 case we have the antisym- 
metric combination 



which gives the weight 5 trigonal PDE 

pnu = 6 pnpu + 3 Hspii- 

For the symmetric case we get a derivative of the weight 4 equation, plus (i 
multiplied by the same equation. With even (odd) weights, the symmetric 
(antisymmetric) combinations give the 4-index pij^i relations. 
At weight 6 we have the three symmetric diagrams 



which give a set of three overdetermined equations with the unique solution 

Pill = ^Pll + Pl2 + 4 Pl3 - 4 P11P22, 

P1122 = 4 pi3 + 2 /ig + 4p^2 + 2 P11P22 + 3 H3pi2- 

Continuing in this way we recover the strictly trigonal versions of the full set 
of equations given in |12] . 

6. Derivation of the integrable hierarchies via residues 

This method of deriving integrable equations on the Jacobians of algebraic 
curves is based on the following relation due to Sato, [28], [29], as applied to 
the algebro-geometric context by Fay, (T5] . 

The Hirota bilinear equations 



Res5=o|^ 



00 



CO 



(6.1) 



n=l J I n=l 

00 ^ f OD 



e^p E<^^" r^p -E 




0, 



. n=l J (_ n=l 

where ^ = ^{Q) is the local coordinate of the point Q near P, ^(P) = defining 
the KP hierarchy are equivalent to the following Bilinear Identity 
Denote ^{f) the differential 



(6.2) 



k=l 
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where fifc(P) are normalized second kind differentials with k-th order poles at 
infinity, so that Q{t, P) has an essential singularity 



and the ^lk{P) are defined by 



(6.3) "=1 

/ nfc(P) = 0, j = l,...,^7. 

Theorem 6.1 (Bilinear Identity). Let be a small circle about the point 
p G X, then 

j^^ 9 (^{Q) - 2l(P) - - ^ j ^ (^^^(Q) - 2l(P) - ^^^^ + ^) 

(6.4) xexpjy fi(t + t',P)|E-2(P,g)de(g) = 0, 

where 2l(P) = /^dv - Abelian image of a point P, ^{Q) is the local coordinate 
of the point Q near P, ^{P) = 0, E{P,Q) is the Schottky-Klein prime form 
Ii2.9\) and the integration is around the unit circle centred at the point p, = 

{g,IC(g)l = i}. 

Proof. Note that the integrand in (16.41) is holomorphic in q in the interior of 
5*^ except at its centre Q = P where it has an isolated essential singularity 
with zero residue. Therefore this integral vanishes. 

To prove the equivalence of (16.11) and (16. 4p note first that 

(6.5) n{t + t') = n{t) + n{t'). 

In ([63D, the factor 



(6.6) 



e (^A{Q) - AiP) - Uiti - n j exp I j Q{t) 



can be expressed as 
(6.7) 

tnC" > d AiQ) - A{P) - Uiti - U 1 exp j Y ^rnntn^ \ • 

n=l ) \ i=l / \rn,n=l 
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Taking into account the form f l3.2p for r(t; u), the last two factors in the above 
formula may be expressed 
(6.8) 

A{Q) - A{P) - UiU - w I exp <^ J2 ^mntn^ ^ = exp <^ - J] ^— }> r(t; n), 

1=1 / \'m,n=l ) K. n=l " 

which completes the proof of the equivalence of (16 ■4p and (16. ip . 

Remark 6.2. Equations of this form were first written down in terms of 'ver- 
tex operators' 

exp|^t„x-|exp|-f;i^A 

for general r-functions in the work of Sato. They may be understood as gen- 
erating functions for integrable hierarchies of Hirota bilinear equations. 

Using this theorem we may obtain partial differential equations relating the 
Kleinian symbols, pij, pijk etc. To do this we now substitute the r-function 
(11.81) into this expression and compute the residue. This is parameterized by 
t and e. The coefficients of monomials in t are differential-difference expres- 
sions in a at arguments -\-u and —u. To get purely differential expressions, 
we now let n — )■ 0; the 'time' derivatives act on a via its n- dependence, 
though, so o"j(w) = — (Tj(— n) — )• ±(Tj(0), etc. We then replace the derivatives 
of the (T-function by derivatives of the p function using the recursive relations 
(13.51) described earlier. This recovers the partial differential equations for the 
Kleinian functions which we need. 

Example I 

In the genus 2 hyperelliptic case, the first nonvanishing relations occur at 
t-weight 3 in the tj, i.e. the coefficients of tf or t^. Both these give the first 
4-index relation 

pun = Qpn + «4pii + 4pi2 + |a3- 
The next non-zero term is at weight 5, where we recover the relation for piu2, 
etc. In contrast with the previous section, this approach only gives the even 
derivative relations at odd t-weights. 

The calculations by this approach, because they involve bilinear products of 
r-functions, soon become very computer-intensive, and we have not pursued 
them very far. 

Example II Within this method we are able to recover all quadratic relations 
for four indexed symbols pijki, ^ > i, j, k,l < 3 and cubic relations for three 
indexed symbols p'^j^- We report here the quartic relation between even vari- 
ables, which must be one of the relations defining the Kummer variety of the 
(3,4)-curve, 

Pl2 - P22 - 2 pup33 + 2pl^ + 4 + P1II3P22 - QpUpl3p22 + ^P\lpYlp2Z 
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+ P11P22 - 2piip^2p22 - |plll3p?i + 8pnpl3 + 2/ii2 + 4;U6P?1 - 4/i6pllp22 
+ 3/i3pl2pl3 + 3/i3pllp23 + /U3P?2 + 2 /LteplS + /"epL + /W9P12 - /U3^P?1 
- 3yU3piipi2p22 = 0. 

This relation is of independent interest since it is of weight 12 and cannot be 
written in the form flS.lOp . The detailed structure of the Kummer variety in the 
(3,4) case is receiving further investigation and will be reported on elsewhere. 

7. Discussion 

In view of the "compare and contrast" objectives of our paper, the reader 
may wish to know which of these methods is the most effective when calculat- 
ing the required PDEs. Unfortunately this is not an easy question to answer. 
Currently the two methods associated with the r-function take rather longer 
to execute, with the residue method in particular being slow and with large 
memory overheads. However this may be due in part to the fact that we have 
extensive experience over ten years or more in developing computer algebra 
code for the "Classical method" in all its variations. We have much less ex- 
perience in working with the tau-function methods. So it may be that with 
more study and with more computational experience, the r-function methods 
become more competitive. We would stress, however, that both r-function 
methods are more systematic than the classical methods, and may provide a 
suitable way of calculating the PDEs associated with more complicated curves. 
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